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Abstract 

The direct production rate of ip in the T decay is shown to be dominated by 
the process T — > ggg* followed by g* — > ip via the color-octet mechanism proposed 
recently to explain the anomalous prompt charmonium production at the Tevatron. 
We show that this plausibly dominant process has a branching ratio compatible with 
the experimental data. Further experimental study in this channel is important to test 
the significance of the color-octet component of cc pair inside the ip system. 



internet address: cheung@utpapa.ph.utexas.edu 
2 Internet address: keung@fnalv.fnal.gov 
3 Internet address: yuantc@ucdliep.ucdavis.edu 



1 Introduction 



The most appealing explanation of the excessive production rates of prompt ■0, and Xcj 
observed at the Fermilab Tevatron [ffl, Q| is given by the combination of the ideas of gluon 
fragmentation into quarkonium || and the color-octet mechanism 0] , in which a gluon frag- 
ments into a color-octet 3 Si cc pair which subsequently evolves nonperturbatively into the 
physical charmonium states by QCD dynamics. While the nonperturbative parameters as- 
sociated with the color-octet mechanism must be extracted phenomenologically from the 
rates of prompt charmonium production, the prediction of the shape of the transverse mo- 
mentum spectrum agrees well with the data 0- 0- ^ 0- A comprehensive review of these 
two theoretical developments and their implications at the Tevatron and LEP can be found 
in Ref. ||. If this mechanism is correct, it may give rise to many testable predictions for 
charmonium production in Z° decay [jj, low energy e + e~ annihilation [TO], photoproduction 
at fixed target and HERA experiments ||11|| , hadroproduction at fixed target experiments 
T2\ and at LHC O], and .B-meson decays [IJ]. Double prompt quarkonium production 



from the color-octet mechanism has also been studied at the Tevatron [15]. In this paper, 
we show that the color-octet mechanism can also provide the dominant contribution to ip 
production in T decay. 

The available experimental data on charmonium production in T decay are listed as 
follows, 

= (1.1 ± 0.4) x 1(T 3 CLEO [16], 

< 1.7 x 1(T 3 Crystal Ball |TJ, (1) 

< 0.68 x 10~ 3 ARGUS RIB], 



Br(T -> ip + X) 



in which the CLEO and the Crystal Ball results show a slight inconsistency. But it is 
interesting to note that these experiments can reach the branching fraction for T — > ip + X 
at the level of 10 -4 — 10~ 3 . Trottier [19| studied the indirect ip production in T decay via the 
production of intermediate physical Xc states, which decay radiatively into ip. However, this 
indirect production of ip contributes to a branching ratio less than 10~ 4 in T decay. Since 
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the experimental branching fraction is already at the level of 10~ 4 — 10~ 3 , we definitely need 
a new production mechanism to explain the data. We also hope that the slight discrepancy 
in the above experiments can be resolved in the near future. 

Conventional wisdom tells us that hadronically T decays predominantly through bb an- 
nihilation into three gluons. The rich gluon content in the final state makes it rather easy 
for the gluon to split into a cc pair in the color-octet 3 Si configuration. If ip can be formed 
from this color-octet configuration at a significant level as predicted by Braaten and Fleming 
0] at the Tevatron [|l], 0] , charmonium states should be abundantly produced in T decay. A 



previous theoretical study [20] of the process T — > ip + X was based on the color-evaporation 



model [21], with which the color-octet mechanism shares some common spirit, but the model 
fails to be systematic. Another qualitative estimate for T — > ip + X can be found in Ref . [^2 



In Section 2, we will briefly review the description of the inclusive decay and production 
of quarkonium based on the nonrelativistic QCD (NRQCD) factorization formalism given 
by Bodwin, Braaten, and Lepage |2"3]| . In Section 3, we will discuss in detail several new 
color-octet processes relevant to ip production in T decay allowed by the general factorization 
formula. In Section 4, we compare the production rates of different processes and discuss 
the energy spectrum of ip in T decay. 



2 NRQCD Factorization Formalism 

The factorization formalism [2tJ for the inclusive decay and production of heavy quarkonium 
allows us to probe the complete quarkonium Fock space in a systematic and consistent 
manner based upon NRQCD. It can be straightforwardly applied to the case of inclusive 
charmonium production from bottomonium decay |Tj|. For the case of T — > ip + X, we have 
the following factorization formula, 

dT(T ^ ^ + = £ df mn {T\O m \T) {Of) , (2) 

m, n 
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where dT mn are the short-distance factors for a bb pair in the state m to decay into a cc pair 
in the state n plus anything, where m, n denote collectively the color, total spin, and orbital 
angular-momentum of the heavy quark pairs. dT mn can be calculated in perturbation theory 
as a series expansion in a s (m c ) and/or a s {m^). Contributions to dT mn that are sensitive to 
the quarkonium scales (m^v b or m c v c or smaller, where v b and v c are the relative velocities of 
the heavy quarks inside the bound states) and to Aq C d can be absorbed into the NRQCD 
matrix elements (T|O m |T) and (0%). We use the notation (0%) to denote the vacuum 
expectation value (0|O^|0) of the operator 0%. The nonperturbative factor (T|O m |T) is 
proportional to the probability for the bb pair to be in the state m inside the physical bound 
state T, while (0%) is proportional to the probability for a point-like cc pair in the state n 
to form the bound state ip. The relative importance of the various terms in the above double 
factorization formula (0) can be determined by the order of v b or v c in the NRQCD matrix 
elements and the order of a s in the short- distance factors dT mn . 

In the color-singlet model [f£3}| , the NRQCD matrix elements involved in the process 
T — ► ip + X are (T|Oi( 3 S'i)|T) and (Of( 3 Si)). According to the velocity scaling rules 
P3| , they are scaled as mfff and m 3 f 3 , respectively, and can be related to the quarkonium 
wavefunctions as follows: 

(TldC^OlT) « §I^t(0)| 2 , (3) 
(OfCS,)} « 3^|^(0)| 2 , (4) 

with N c denotes the number of colors. Therefore, these color-singlet matrix elements can 
be determined from the leptonic widths of the T and ip. The short-distance factor in the 
color-singlet model for this direct process includes bb( 3 Si, 1) — > cc( 3 Si,l)gg and bb( 3 Sx, 1) — > 
cc( 3 Si,l)gggg. The two possible color configurations of the heavy quark pair are denoted by 
1 for singlet and 8 for octet. The leading order Feynman diagrams for these processes are of 
order a®. Due to such a high order in the strong coupling constant, it is unlikely that these 
color-singlet processes can be the dominant production mechanism. 
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The first example of the double factorization formula such as (0) is the indirect ip pro- 



duction in the decay T — > Xcj + X with Xcj ^7 considered by Trottier [Hi]. The 
short-distance factor r(&6( 3 5i, 1) — > cc( 3 Pj, 1) + ggg) is of order c^, which is enhanced by 
a factor of l/a s compared with the direct color-singlet processes mentioned in the previ- 
ous paragraph. However, the infrared divergence in the leading order calculation of the 
short-distance factor indicates that the results are sensitive to the scale m c v c or smaller. 
Therefore, in addition to the color-singlet matrix element (0* C,J ( 3 P/)) (scales like m c vi), one 
also needs to include the color-octet matrix element (0g CJ ( 3 Si)) (scales like m 3 t>f) to absorb 
the infrared divergence. The short- distance process associated with the color-octet matrix 
element is bb( 3 Sx,V) — > cc( 3 Si, 8) + gg, which is of order aj. In this case, the introduction 
of the color-octet matrix element is required by perturbative consistency, since the infrared 
divergence would otherwise spoil the one-term factorization formula. 

In the next section, we will consider the direct and indirect ip production in T decay 
with short-distance factors of order a 3 and aj.. These are possible only if higher Fock states 
of the ip or T are considered. We will also consider processes with short-distance factors 
of order aa 3 , a 2 a s , and a 2 a 2 s that are suppressed by electromagnetic coupling but may or 
may not require higher Fock states of the quarkonia. In the following, we will first consider 
the produced cc pair in the color-octet 1 5'o, 3 S\, or 3 Pj configuration, which subsequently 
evolves into physical ip described by the matrix elements {Of ( 1 So)), (Of ( 3 Si)), or (Og( 3 Pj)), 
respectively. These color-octet matrix elements are suppressed by v\ relative to the color- 
singlet matrix element (Of ( 3 Si)). The matrix element (Og ( 3 Si)) has been extracted from 
the CDF data 0, [5], p|, while two different combinations of the other color-octet matrix 
elements have been extracted from the CDF data M and from the photoproduction data by 
Amundson et al [IT|. 

Though of much smaller effects, we also consider the contributions by the higher Fock 
state of the color-octet bb pair inside the T associated with the matrix element (Y|08( 3 5i) |T), 
whose value has not yet been determined. An order of magnitude of this matrix element 
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can, in principle, be estimated by considering the ratio 

(T|Q 8 ( 3 g 1 )|T) (gfgft)) (vj 
(T\0 1 ^S 1 )\T) (Of^S,)) ~ \vl 
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(5) 



which implies that its value should be highly suppressed. The ratio in (||) tells us that 
processes associated with a color-octet cc pair inside the produced ip and a color-singlet bb 
pair inside the decaying T are much more important than those with a color-octet bb pair 
inside the T and a color-singlet cc pair inside the ip. Using the value (Ot{ 3 Si)) « 0.014 
GeV 3 |, g |, (Of( 3 Si)) « 3(^|Oi( 3 Si)|^) « 0.73 GeV 3 from the leptonic width of *p, 
(T|Oi( 3 5i)|T) « 2.3 GeV 3 from the leptonic width of T, and v 2 c « 0.3 and v 2 « 0.08, we 
obtain (Y|08( 3 Si)|Y) « 3 x 10~ 3 GeV 3 . However, such a large value for this matrix element 
would substantially increase the hadronic width of T, which would diminish the leptonic 
branching ratio to an unacceptable level. Obviously, this matrix element enters into the 
hadronic width of the T via the short-distance process 66 ( 3 Si, 8) —>■ g* —>■ qq. In order not 
to spoil the experimental value for the leptonic branching ratio and the total hadronic width 
of T, it is necessary to put a bound on the value of the matrix element (T|Og( 3 Si)|T). The 
hadronic width of the T has the factored form 



r(T -> light hadrons) = [ f (66( 3 Si, 1) ggg) + £ f(66( 3 Si,l) 

xOTIO^SOIT) 



q=u,d,s,c 

/3, 



+ £ f(66( 3 Si,8)^gg)(T|0 8 ( 3 Si)|T) + --- , (6) 

q=u,d,s,c 

with the following short- distance factors calculated to leading order in a and a s , 



20a 3 
h 



r(66( 3 S 1; 1) - ggg) = (vr 2 - 9) , (7) 



r(66( 3 Si,i)^//) - - 

and 



27iN c QlQ 2 f a 2 
3m? 



7Ttt 2 



T(bb(- i S 1 ,8)^qq) = ^ . (9) 



'b 
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In Eq. (§), N c is 1 and 3 for / equals charged lepton I and light quark q, respectively; Qf 
is the electric charge of the fermion / in unit of the positron charge; and we have set m q 
and mi to zero in Eqs. (||) and (Qj for simplicity. Using the following expression for muonic 
branching ratio 

BR(T - ) = r(T - light hadronsf+ EJtf FFj ' 



together with the experimental value for BR(T — > ) = 0.0248 ± 0.0007 |25 

obtain the following bound on (Y|08( 3 Si)|T): 



we can 



(T|0„( 3 S,)|T) 



( 

1.9 



x 10" 4 GeV 3 , (11) 



where we have allowed a 2a variation on BR(T — > fi + {i~). Alternatively, we can obtain 
another bound by using the total width of T, but the result is not as good as the one given 
by Eq. (|Tl"D . Therefore, in the rest of the paper we will use the value (Y|0 8 ( 3 Si)|Y) ~ 5xl0~ 4 
GeV 3 . With this value for (Y|08( 3 Si)|Y), we obtain the ratio 

f Wft.g) ^g*^qg)(r\o 8 (3s 1 )\r) Qoi7 

fibb^l) ^ gggXTlO^lT) " ' ' 
which is now sufficiently small posing no threat to the experimental value of the leptonic 
width in the Y decay. The input parameters used in our later numerical calculations are 
summarized in Table [I| for convenience. Since all our new calculations are at tree-level only, 
we, to be consistent, extract the color-singlet matrix elements from the leptonic widths of 
if), if)', and Y using tree-level formulas. 



3 Color-Octet Processes 

We shall first consider two processes with a color-octet cc pair inside the produced ip and a 
color-singlet bb pair in the decaying Y. 
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3.1 66( 3 5i,l) 



r 



cc{^Lj,8)+g 



The leading order diagrams for the process bb( 3 Si, 1) — » cc( 2S+1 Lj, 8) + (7 is of order a 2 a s , 
one of which is shown in Fig. |l](a). This is similar to the process Z —>■ cc( 2S+1 Lj,8) + g, 
which is negligible because the short-distance factor is suppressed by powers of m 2 /M| from 
the quark propagator ||. But in the present case it is only suppressed by powers of m 2 /m 2 . 
We will restrict ourselves to the case of L = and 1 only. Although the contributions from 
higher values of L can be included easily, they are further suppressed by powers of v 2 . The 
inclusive production rate from these processes can be written as 



r la (T^ + X) = (TlOifSOIT) 



(O 8 YS )>WSi,I) - ccCS ,8)+g) 



+ E (Oi( 3 Pj)}T(bbCS u l)^ccCPj,8)+g) 

J=0,l,2 

The short-distance factors are calculated to leading order and are given by 



r(fofo(^i,I)^cc( 1 i S'o,8)+^) 
f(66( 3 S 1 ,l)^cc( 3 P ,8)+^) 
f{bb{ 3 S 1 ,l)^cc( 3 P 1 ,8)+g) 
f(bbCS 1 ,l)^cc( 3 P 2 ,8)+g) 



4n 2 Q 2 c Q 2 a 2 a 



m\m c 



;i-o, 

m^m 3 1 — £ 

1 i + e 



Art 2 Q 2 c Q 2 a 2 a s 
9 

87T 2 Q 2 c Q 2 b a 2 a s 

9 mfm 3 1 — £ ' 

8ix 2 Q 2 c Q 2 a 2 a s 1 1 + 3£ + 6£ 2 
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(12) 

(13) 
(14) 
(15) 
(16) 



where £ = M|/M| « m 2 c /m 2 . We have also used the nonrelativistic approximation for 
the mass of the quarkonium: Mx ~ 2m& and « 2m c . We note that &6( 3 5i,l_) — > 
7* — > cc( 3 S'i, 8) + g vanishes. Using the heavy quark spin symmetry relation (Og( 3 Pj)) ~ 



(2 J + l)(Og ( 3 Pq)) (HI) the total width from these processes can be simplified as 



r la (r^i> + x) 



An 2 Q 2 c Q 2 a 2 a s (TIO^OIT) 



3m 2 



m\m c 



(otCs ))(i-o 



(1-30 2 6(1 + 2(l + 3£ + 6£ 2 



1-e 



+ 



1-e 



(17) 



which can be normalized by the leptonic width of T: 

2nQ 2 a 2 (T\0 1 ( 3 S 1 )\T) 



r(T 



e e 



m 



(1* 



Using m c = 1.5 GeV, m b = 4.9 GeV, a s {2m b ) = 0.179, (C^(%)) « (C^( 3 P ))/™c 
10" 2 GeV 3 0, and BR(T -> e + e~) « 0.0252 @, the contribution from the above color- 
octet processes to the inclusive branching ratio BR(T — > ■?/> + X) is only 1.6 x 10~ 5 , which is 
almost two orders of magnitude below the CLEO data (|l]). 

3.2 bb{ 3 S h l) ggg* cc( 3 S 1} &) + gg 

Fig. |T](b) shows one of the six Feynman diagrams for the bb( 3 Si,V) pair annihilating into 

three gluons with one of the gluons converting into the cc( 3 Si,8) pair. This process is of 

order a A s and its calculation is very much similar to the process T — > gg^y* — > gg 11 ||26|| . 

Introducing the following scaling variables: 

Ei E E 

mb m-b mb 

such that x v + xi + X2 = 2, where Ei stands for the energy of the particle i. In the rest frame 

of T, the differential decay width is given by 

dTlb 'T - i>+X) = -^—(bbCSx,!) - cc( 3 5i,a)+^)(T|0 1 ( 3 5i)|T>(Of( 3 5 1 )> , (20) 



dx v dxi dx v dx\ 
with 



x 



^_(66(3 5l , 1) - cc( 3 S u 8) +gg) - J™^ ^ _ ^ ^ _ - _ ^ 



2^ + 2^(6 - Ax v + 2xi - x v X! - x() 
+ 2£ 2 (ll - 16z„ + 6x 2 v - (8 - 2x c - x 2 )x! + (4 + x v )xj) 
+ £{4(1 - x v ){A - 5x v + 2a: 2 ,) - (32 - 44a:,, + Ux 2 v )x 1 + (20 - 18x v + xl)x\ 

2^2 Xyj'jjC ^ - f~ x 

+ 2 ^2 - 6a;,, + 7a: 2 - 4a: 3 + x 4 v - (6 - 13z„ + 9x 2 - 2x 3 )xi + (7 - 9a:„ + 3a: 2 )a; 2 

2^2 X^^X-^ — r~ X 

where the ranges of integration for x v and X\ are 

lJl< x v <1 + C, (22) 



(21) 



l(2-x v -y/xl-4£) < x, <l(2-x v + ,Jxl-4£) . (23) 



One can integrate over x\ to obtain the energy distribution of ijj, 



dx v K K v '" y yay 486m 3 m 2 (2 - x v f(x v - 2£) 2 

4(8 + 8£ - 14£ 2 - 2£ 3 - \2x v - A£x v + lOf 2 ^ + 5a; 2 - £x 2 - £ 2 x 



,'2-x v - Jxl -4f\ 

x (l + e-x„)log| V , = ) + (2-x,)v/x 2 -4e 

2 - x„ + Jx 2 , _ 4 £ 



x (16 + 28£ + 20£ 2 + 4£ 3 - 24x, - 36fx w - 12£ 2 x„ + 14x 2 + 13£x 2 - 4x 3 ) . (24) 
Numerically integrating x v , we obtain the partial width, 

r lb (T ^ + X) = ^ {T\0 1 ( 3 S 1 )\T)(Ot( 3 S 1 )) x (0.571) . (25) 



We note that the color-octet piece for the process T — > % c j + X considered by Trottier |19| 
can be obtained from the above Eq. (j25l) by simply replacing the matrix element (Og ( 3 £i)) 
with (Og CJ ( 3 S'i)). However, the expressions for the energy distributions (21|) and ([24]) were 



not given explicitly in Ref. [ I9| . The prediction of the partial width is sensitive to the 
values of the two NRQCD matrix elements, the running coupling constant a s , and the heavy 
quark masses. For convenience we can normalize this partial width to the three-gluon width 
r(T — > ggg) given by Eq. ([7]) times the matrix element (Y|Oi( 3 Si)|Y): 

T lb (T^ij + X) na s (Ot( 3 S 1 )) 0.571 



T(T — > ggg) 8 m 3 n 2 — 9 

It is insightful to take the scaling limit of m b — > oo with x v = E^,/m b held fixed in Eq.(|2 
In this scaling limit, one can deduce 

T lb (T^ij + X) na s (Ot( 3 S 1 )) 



(26) 



(27) 



r(T -> 8 m 3 

The right-handed side of Eq. (p7|) can be recognized to be three times the gluon fragmen- 
tation probability P g -+^ in the color-octet mechanism obtained by Braaten and Fleming [Q. 
Thus, the above scaling limit corresponds to the fragmentation approximation. Comparing 
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the above limit with the exact result (0), we see that fragmentation is not a good 
approximation. However, this limit suggests that the scale to evaluate the a s in Eq. ( p6| ) 
should be 2m c instead of 2m b . With a s (2m c ) = 0.253, (0^( 3 Si)) = 0.014 GeV 3 @, §, and 
assuming BR(T — > (7(7(7) « BR(T — > light hadrons) = 0.92 f2"5| , we obtain 

r lft (T->y> + x) r lfc (r^ + x) 

7^7 = — x BR T -> ^ 

r to tai(T) r(T ->• 

« 2.5 x 10~ 4 . (28) 

This prediction is smaller than the CLEO data by merely a factor of 4, and is consistent 
with the bounds from Crystal Ball and ARGUS. 

The color-octet process studied in this subsection applies to the case of if)' as well, sim- 
ply by replacing the matrix element (Of ( 3 S\)) with (0$ ( 3 Si)), whose value has also been 
extracted from the CDF data to be 0.0042 GeV 3 g §, §, [|. With BR(ip' -> if) + 7) « 57% 



25| , we obtain a branching ratio of 4.3 x 10~° for the inclusive ip production in the T decay 
that comes indirectly from ip'. 

One can also consider the processes T — > 7(7(7* followed by g* — > if) (ip r ) via the color-octet 
mechanism, and T — > gg^y* followed by 7* — > if) {if)') in the color-singlet model. Up to overall 
normalization, the energy spectrum of the if) for these two processes are predicted to be 
same as in Eq. fl2"4j) . However, their partial widths are suppressed by factors of 8a/(15a s ) ~ 
0.02 and 32a 2 (Of( 3 S 1 ))/(A5a 2 s (Ot( 3 S 1 ))) ~ 0.06, respectively, compared with the width of 
Eq. (P5]) . Thus they contribute a branching fraction about 2 x 10 -5 in the inclusive decay 
T — > ip+X. The indirect contribution from the if' from these two processes is about 6 x 10~ 6 . 



3.3 66( 3 5' 1 ,8)^/^cc( 3 P J ,l) + r/ 

We now turn to the case where the bb pair inside the T is in a color-octet 3 S\ state. The 
first process we consider is the production of Xcj from T decay. The leading order Feynman 
diagram is depicted in Fig. [l|(c). The factorization formula for the decay rate can be written 
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as 



ric(T - Xcj + X) = r(fofo(^i,8) - cc( 3 P J; 1) +p)(T|0 8 ( d S 1 )|T)(0^( d P / )) . (29) 

We note that up to coupling constants, color factors, and NRQCD matrix elements, these 
processes are similar to the ones bb( 3 Si, 1) — > 7* — > cc( 3 Pj, 8) + g considered in Section 3.1. 
The short-distance factors can be extracted easily from the previous calculations, 

f(bbCS u 8)^cc( 3 P ,l)+g) = ^-^^7^, (30) 



r(66( 3 S 1 ,8)-cc( 3 P 1 ,l)+^) = ci - ^ i 3 ' ( 31 ) 

81 mim'l 1 — 4 

f(66( 3 5 1 ,8)-cc( 3 P 2 ,l)+^) = ^ * 1+ 5 + 6 * 2 . (32) 



51 m 4 m 3 1 — £ 
1 

51 m 4 m 3 1 — £ 
2 a 3 1 1 + 3 
415 m\m\ 1 — £ 

The matrix element (0i cJ ( 3 .Pj)) is related to the wavefunction according to |23] 

, , , Nl2 2vr (Of J ( 3 Pj)) , . 

p; 2 « — / J;/ . 33 

9 2,7+1 v ; 

Using the value of the matrix element (T|08( 3 5'i)|T) « 5 x 10 -4 GeV 3 estimated in the last 
Section, a s (2m b ) = 0.179, |P^(0)| 2 = 0.075GeV 5 from the potential model calculation [|7|, 
and the branching ratios BR(x cl — ► ip + 7) = 0.273 and BR(x c2 — > ip + 7) = 0.135 (xco has 



a negligible branching ratio into ip) [p5 |, we obtain the width T lc (T — > ?/> + X) « 0.05 eV 



and thus a branching ratio of 9 x 10 7 . Therefore, these indirect contributions are negligible 



when compared with the indirect mechanism considered earlier by Trottier [19 . 



3.4 66( 3 Si,8) g* cc( 3 5i, 1) +^ 

One of the six leading Feynman diagrams for the short- distance process bb( 3 S\,8) — > g* — > 
cc( 3 S'i, 1) + (7(7 is shown in Fig. |l](d). The corresponding differential width can be written as 
the following factored form 

dVld (T — > *0 + Jf) = -^—(tof^gj - cc( 3 S 1 ,l)^)(T|0 8 ( 3 5 1 )|T)(Of( 3 S 1 )) (34) 



doc d dec \ dec u dec i 
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with the short- distance factor given by 

<,l: -(66(^,8) -+cc( 3 S 1 ,l)<?<?) '^^ 



dx v dx\ ' 486m^m c (2 — x 1) ) 2 (l — ^ — Xi) 2 (l + ^ — Xt, — Xi) 2 

x + 2£ 3 (8 - 5x v + x 2 v - 2xi + x v X! + x\) 

+£ 2 (^28 - 46x„ + 2lx 2 v ~ &4 - (28 - 26x v + 6x 2 )xi + (14 - 6x v )xj 
+2^6 - 17x v + I6x 2 v ~ 6^ + xi~ (12 - 22x,, + 12X 2 - 2x 3 v ) Xl 
+ (10 - 12x v + 3x 2 )x 2 - 2(2 - x v )x\ + x 4 ^) 

— (1 — xi)(l — x v — xi){l — x v + 2x\ — x v X\ — x 2 )j , (35) 



where the allowable ranges of x v and X\ are given in Eqs. (|22|) and (|23|) . Integrating over x 1; 
we obtain the energy distribution for ip 

dV . l7/ o„ oN /Sr , lN . 57ra 4 1 1 
-(&&( 3 ^i,S) -*cc( 3 S u l)+gg) - 



dx v ' 1 " '' 486 mfm c (2 — x v ) 2 (x v — 2£) 3 



x 



| (4 + 20£ + 28£ 2 + 16£ 3 - 12x v - 36£x,, - 24£ 2 x J3x 2 + 14£x 2 - 4x 



2i-x v + \X 2 V -Ai 



x(x„-20 ^/x2-4^-41og , 

2^-^-^x2-4^ 
(2f + 16£ 2 + 6£ 3 - 16£ 4 - 8£ 5 - 12£x„ - 20£ 2 x„ + 24£ 3 x 1 , + x 2 

+20£ 4 x, + 12£x 2 - 8£ 2 x 2 - 17£ 3 x 2 - x 3 - £x 3 + 5£ 2 x 3 ) } . (36) 



The partial width is obtained numerically: 

r li( r + = ^ « M£M x (1 , 3 o, . (37) 

486 mi m c 

With the previous inputs, a s (2mb) = 0.179, and (Of( 3 Si)) ~ 0.73 GeV 3 obtained from the 
leptonic width of ip, we obtain a width of 0.02 eV only, which gives a branching ratio of 
3 x 10- 7 . 



3.5 Other Color-Octet Processes 

When T annihilates into a qq pair via the s-channel photon 7*, a bremsstrahlung virtual 
gluon emitted from the light quark line can become an octet cc( 3 Si,8), which then turns 
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into if). The factored form of the rate for this process is 

T le (T - if; + X) = f(66( 3 5 1; l) - cc( 3 S 1 ,8) + qq)(T\0 1 ( 3 S 1 )\T) (Ot( 3 S 1 )) . (38) 

This process is of order a 2 a 2 sl which is similar to a potentially dominant color-octet process 
in the prompt if) production in Z° decay [BJ. We can easily translate the formula from Ref.[H 
to the present case. The partial width is given by 

Tle(T - 7 * - <?#) a ^ 2mc )^(^l))f 5(1 _ e2) _ 2eloge , \ nTi ( Z 



T(T — > 7* — > gg) 36 m 



2Li 2 



1 + t 



- 2Li 2 f ^—-^J -21og(l + 01oge + 31oge + log 2 e 







(39) 



Here Li 2 (x) = — log(l — t) is the Spence function and £ = [m c jm\ } ) 2 . T(T — ► gg) is 
given by Eq.(j|) times the matrix element (Y|Oi( 3 Si)|Y). Using mj = 4.9 GeV and m c = 1.5 
GeV, £ ~ 0.0937 and the curly bracket in (p9j) is about 1.1, as compared to the much larger 



value of 27.9 in the corresponding case in Z° decay || where £ = m^/M| w 1.1 x 10 3 . 



Therefore, the double logarithmic terms in Eq.(p9|) do not provide sufficient enhancement in 
our present case. Numerically, the ratio in (|3^) is only 8.3 x 10 -6 , and so this process can 
be safely ignored. 

In addition to the above processes, one can also consider a color-octet 3 Si bb pair an- 
nihilating into a light quark pair via a s-channel gluon, followed by an off-shell photon 
bremsstrahlung off either the b quark line or the light quark line and eventually turns into 
the if). This process also involves the matrix element (T|Og( 3 S'i)|T) and, therefore, will 
be suppressed. One can also consider both the bb and cc pairs inside the quarkonia are in 
the color-octet configuration. Processes associated with such configuration are further sup- 
pressed by powers of v 2 and v\ compared to those we have studied in this paper. We ignore 
them in this work. 
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4 Discussions and Conclusions 



The energy distribution of charmonium in the T decay can provide an interesting test for 
the NRQCD factorization formalism discussed in Section 2. The energy distributions for 
the processes of Figs. 0(a) and 0(c) are just a delta function with the peak at one half of 
the T mass, while the energy distributions for the processes of Figs. |l|(b) and |l](d) are given 
by Eqs. (|2~3| ) and PB]), respectively, and are shown in Fig. ^. The solid curve in Fig. ^ is 
the energy spectrum of if) for the dominant process &&( 3 5i,l) — > g*gg — > cc( 3 Si, 8)gg; it is 
monotonically increasing as the if) energy increases and eventually cut off by the kinematic 
limit. In the fragmentation approximation (m& — ► oo with E^/rrib held fixed), the energy 
distribution of if) for this subprocess is the same as that of the fragmenting gluon, up to an 
overall constant. The result using the fragmentation approximation is shown as the dashed 
curve in Fig. for comparison. The fragmentation approximation has overestimated the 
exact result for all energies of the if). On the other hand, CLEO obtained a relatively 
flat momentum spectrum. However, the rising spectrum shown by the solid curve in Fig. |2] 
can be softened by various mechanisms such as final state interactions of the soft gluons 
P%| and by relativistic corrections of the bound state PDfl , just to name a few possibilities. 
This situation is very similar to the photon spectrum in the decay T — ► jgg pUf . 

The dashed-dotted curve in Fig. || is the energy distribution for the suppressed channel 
bb( 3 S%,8) — > g* — > cc( 3 Si, l)gg and it has a much flatter shape. 

Among the several processes that we have considered in this paper, those two with the 
color-octet bb pair inside the decaying T are more than two orders of magnitude below the 
CLEO data and therefore negligible, while the other with the color-octet cc pair inside the 
produced if) are relatively more important. The largest contribution is being the process 
shown in Fig. |l|(b), which has a branching ratio about 2.5 x 10~ 4 . The next largest contribu- 
tion is the indirect process considered by Trottier From Eq.(21) of Ref. fl9|] and using 



the heavy quark spin symmetry relation (0g cj ( 3 5'i)) ~ (2J+ l)(O8 c0 ( 3 S'i)) |23| and the value 



of the matrix element {Of °( 3 Si)) = 0.0076 GeV 3 obtained by fitting the Tevatron data H g, 
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we obtain X)jBRTrottier(Y ~~ ¥ XcJ + X] Xcj —> ~ 5.7 x 1CT 5 . Processes that are compara- 
ble to this one are (1) the process shown in Fig. |l](a) which has a branching ratio of 1.6 x 10~ 5 , 
(2) the processes T — > gg^* — > ip + X and T — > ^gg* — > ip + X via the color-octet mechanism 
which have a combined branching ratio about 2 x 10~ 5 , and finally, (3) the indirect contribu- 
tion from ip' having a branching ratio about 5 x 10~ 5 . Therefore, adding up the contributions 
from all these processes, we obtain a branching ratio BR(T — > ip + X) « 4 x 10~ 4 , which is 
within 2cr of the CLEO data. Given these theoretical results, it would be very interesting 
to have more precise measurements of the inclusive rates and energy spectra of charmonium 
from the T decay. This would provide a crucial test of the NRQCD factorization formalism 
applied simultaneously to both the bottomonium and charmonium system. 

This work was supported in part by the United States Department of Energy under Grant 
Numbers DE-FG02-84ER40173, DE-FG03-93ER40757, and DE-FG03-91ER40674. 
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Table 1: Input parameters used in our calculations. 



NRQCD matrix elements 


Value 


(Of( 3 5 1 ))^3(^|0 1 ( 3 5 1 )|^) 


0.73 GeV 3 


(Of'eS 1 ))^3(4>'\Oi( 3 SiM') 


0.42 GeV 3 


l^ c (0)| 2 ^f(Of°( 3 P )) 


0.075 GeV 5 |27| 


(Tld^OIT) 


2.3 GeV 3 [Eq.©] 


(ot( 3 s x )) 


0.014 GeV 3 0, |, § 




0.0042 GeV 3 § |, § 


(OtCSo)) « (O 8 ^( 3 P ))/m c 2 


10- 2 GeV 3 § 


(or( 3 50) 


0.0076 GeV 3 § 


(TlOg^OIT) 


5 x 10" 4 GeV 3 [Eq.flll])] 


Other parameters 


Value 


m c 


1.5 GeV 


m b 


4.9 GeV 


a s (2m c ) 


0.253 


a s (2m b ) 


0.179 



19 



Figure Captions 



1. Some of the contributing Feynman diagrams for the short-distance processes: 

(a) bb( 3 S 1 ,l)^ 1 *^cc( 2S + 1 L J ,8)g ] 

(b) bb( 3 S 1 ,l)^g*gg^cc( 3 S 1 ,8)gg; 

(c) bb^SuS) ^g*^ cc( 3 Pj,l)g; and 

(d) bbfS^^g^cefSuDgg. 

2. The inclusive energy spectra dT /dx v of ip in the decay T — > ip + X with the followin. 
color-octet processes : bb( 3 Si, 1) — ► — > cc( 3 S\,8)gg (solid) as predicted by Eq. 
compared with the fragmentation approximation (dashed), and bb( 3 S\,8) — > g 
cc( 3 Si,l)gg as predicted by Eq. fl3D|) multiplied by 300 (dashed-dotted). 
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